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Abstract — Wireless sensor networks (WSN), i.e. networks of 
autonomous, wireless sensing nodes spatially deployed over a 
geographical area, are often faced with acquisition of spatially 
sparse fields. In this paper, we present a novel bandwidth/energy 
efficient CS scheme for acquisition of spatially sparse fields in 
a WSN. The paper contribution is twofold. Firstly, we intro- 
duce a sparse, structured CS matrix and we analytically show 
that it allows accurate reconstruction of bidimensional spatially 
sparse signals, such as those occurring in several surveillance 
application. Secondly, we analytically evaluate the energy and 
bandwidth consumption of our CS scheme when it is applied to 
data acquisition in a WSN. Numerical results demonstrate that 
our CS scheme achieves significant energy and bandwidth savings 
wrt state-of-the-art approaches when employed for sensing a 
spatially sparse field by means of a WSN. 

\ Index Terms — Compressive sampUng, sensor network. 

I. Introduction 

Wireless sensor networks (WSN) consist of autonomous, 
, cooperative sensors spatially deployed over a geographical 
area, with applications ranging from surveillance [1] and 
.localization systems [3], [4], to environmental monitoring for 
'physical field sensing and disaster prevention [5], [6]. WSN 
' nodes typically acquire the data and communicate them to a 
node named Fusion Center (FC), which stores the sensors' 
readings or forwards them through wired network infrastruc - 
■ tures for further processing. 

] The availability of energy efficient algorithms for data 
gathering towards the FC is particularly relevant when the 
monitoring network is deployed on a large geographical area 
(e.g. a forest), where highly efficient routing protocols are 
required for sustainable network life-time. Energy efficiency is 
also relevant in those environments where battery recharge or 
substitution may be unworkable (e.g. in underwater networks) 
[7], [8]. On the other hand, efficient exploitation of the avail- 
able bandwidth is an important concern in bandwidth-limited 
or interference-limited environments. 

The unifying sampling and compression approach of Com- 
pressive Sampling (CS) [10] is definitely well-suited to re- 
source limited WSNs' applications. CS based techniques for 
energy efficient WSN data gathering have been recently in- 
vestigated, with particular reference to the trade-off between 
the reconstruction accuracy and the data gathering cost [18]. 
A highly efficient approach is provided by Random Sensing 
(RS), where at each observation time only a randomly drawn 
subset of sensors acquires data and transmits them to the 
FC, typically using single hop links. From a theoretical point 
of view, RS and CS share conditions and procedures for 
signal reconstruction. An energy and bandwidth efficient RS 
procedure appears in [21]. 



WSN monitoring applications are often faced with acquisi- 
tion of spatially sparse signals. A typical example is that of 
temperature monitoring sensor networks for anomalous event 
(e.g. fire) detection: in the early stages of abnormal system 
behavior, in which the event is hopefully detected, the field 
is characterized by one or more small spots at levels largely 
different from the surroundings, and it can be modeled as a 
spatially sparse signal. RS schemes, such as those analyzed in 
[21], poorly perform in sampling signals that are naturally 
sparse in the spatial domain, since the actual number of 
measurements required to reconstruct the field increases and 
the RS bandwidth/energy efficiency deteriorates. 

This paper successfully addresses the efficient compressive 
sampling of spatially sparse signals in a WSN. We introduce a 
novel CS scheme that can be realized in a WSN by distributed 
and parallel data gathering schemes, with restrained energy 
and bandwidth consumption for inter-sensor signaling. Specif- 
ically, the main contributions of this work are the followings: 

• we introduce a novel CS matrix, and we analytically 
demonstrate that it satisfies the CS conditions for sparse 
signal reconstruction; 

• we analytically evaluate the performance, in terms of en- 
ergy and bandwidth efficiency, of a WSN data gathering 
scheme based on the herein presented CS matrix; 

« we show that, on spatially sparse fields, our CS scheme 
outperforms selected RS and CS state-of-the-art ones in 
terms of both energy and bandwidth efficiency. 

The remainder of the paper is organized as follows: Sect. II 
describes the WSN system model. Sect. Ill briefs the CS basics, 
and Sect.IV discusses related works on CS/RS acquisition in 
WSN. In Sect.V we illustrate our original CS scheme and in 
Sect. VI we apply it in a WSN and we analyze the related en- 
ergy consumption and bandwidth occupancy. Sect.VII presents 
numerical results assessing that our CS scheme outperforms 
of state-of-the-art approaches in terms of bandwidth/energy 
efficiency. Sect. VIII concludes the paper 

II. System Model and Network Scenario 

Let us consider a physical field represented by a bi- 
dimensional time varying signal s{x,y;t), {x,y) S TZ^ 
monitored through a a grid of N — Ni x N2 sensors deployed 
over a bi-dimensional covered area, being Ni and N2 the 
number of sensors distributed in the horizontal and vertical 
dimension, respectively. A selected sensor collects the data 
from the others and plays the role of FC. An example of such 
a WSN scenario is depicted in Fig.l: the FC is placed at the 
center of the network, and each sensor is placed at distance 
dk,k=l,...,N -1 from the FC. 
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Fig. 1. Sensor network scenario. 



The sensors periodically measure s{x,y;t) and transmit 
their readings to the FC for monitoring. The sensing period 
St is selected to be almost equal to the coherence time Tc of 
s{x,y;t)^. At time tk = to + 5t, the sensor at the location 
{n\d,n2d) acquires the noisy measurement: 



z [ni,n2] = s {nid,n2d;tk) 



(1) 



where, for simplicity sake, we have dropped the temporal 
variable and we have set the horizontal and vertical inter- 
sensor distances equal to d. Then, the sensors implement a 
suitable dissemination protocol to forward the measured value 
to the FC. The FC collects all the data from the sensors so as 
to reconstruct a representation of the overall field z[ni,n2]- 
Efficient data dissemination towards the FC is widely de- 
bated in the literature, since energy efficiency affects the 
network lifetime, especially relevant in scenarios where the 
deployment of the sensors is difficult or expensive, whereas 
bandwidth efficiency enables WSN monitoring in bandwidth- 
hmited media, e.g. underwater, or in geographical areas where 
different sensor networks coexist. CS and RS paradigms 
provide a theoretical framework for highly efficient field 
monitoring, provided that the monitored data are sparse in a 
suitable domain. We briefly recaU in Secs.III, IV the basics of 
CS and the related work on CS and RS application in WSN, 
respectively. 

III. An Introduction to Compressive Sampling 

Let us compactly represent the bidimensional sequence 
z [ni,n2] by the lexicographically ordered vector z: 

z= [^[l,l]---z[l,iV2] z[2,l]---z[N,,N2]r 

The vector z s is sparse if the number of its non- 
zero samples is restrained compared to its own dimension N; 
rigorously, z is said to be K-spaise if the number of non-zero 
components is K, either in the spatial or in a transformed 
domain (e.g. Fourier, wavelet, etc.). Compressive sampling 
(CS) provides a framework for sensing and compression of 
a sparse signal. 

According to the CS paradigm, compression of sparse 
signals is performed jointly with the acquisition. Specifically, 
z is represented by M projections defined as follows: 



y = *z 



(2) 



being $ a suitable measurement matrix of size M x N . 

A fundamental outcome of CS is that, under suitable con- 
ditions on the sensing matrix if z is ii'-sparse it can be 
accurately recovered from the projections in y provided that 
K < M < N. Specifically, the sensing matrix # must satisfy 
the so caUed Restricted Isometry Property (RIP), i.e., given a 
constant 6, for aU ii'-sparse signals z it must stand: 



(1 



N|^<||*z||^<(l + ^)i|z||^ 



It is proved [13] that, for values of 5 small enough, sparse 
signals can be perfectly recovered from compressive sensing 
measurements. In [30], the authors show that the RIP property 
is verified when the measurements energy ||$i{z|p is more 
and more concentrated, in probability, around the value ||z||2 
as far as the number of measurements increases. 

Reconstruction can be achieved either by solving the fol- 
lowing optimization problem 



argniin ||t||i 



s.t. 



*z 



(3) 



'The coherence time Tc is defined as the time interval over which the 
process almost de-correlates in time. 



or by a greedy iterative pursuit of the support of z; examples of 
this latter approach are provided by the orthogonal matching 
pursuit (OMP) [14] and the compressive sampling matching 
(CoSaMP) algorithms [15]. 

In (2) we recognize that randomly samphng the components 
of z and collecting them in y is equivalent to reaUze CS 
using a particular sensing matrix; therefore, many CS theory 
results apply to RS, too. Both RS and CS techniques have 
been considered for application in WSNs. 

IV. Related Works 

Efficient sensing and data gathering in a WSNs by means 
of CS and RS based techniques has aroused lively interest 
in recent literature. In [11], a CS based distributed, commu- 
nication architecture is exploited to minimize the latency for 
information retrieval under a favorable power-distortion trade- 
off whereas in [20] a CS based sensing and data gathering 
procedure is analyzed for the case of network routing tree 
topologies. In [17], maximization of large scale WSN lifetime 
is pursued by means of a fully distributed algorithm according 
to which each sensor autonomously performs classical or com- 
pressive samphng in order to reduce the number of transmitted 
packets. 

In [18] the authors analyze a RS and multi-hop data gath- 
ering scheme. In this scheme only randomly selected nodes 
measure the field and transmit their readings to the FC through 
the specific multi-hop paths. While a sensor reading is routed 
towards the FC, its value is combined with the ones sensed by 
the sensors in the path, so that each random projection pro- 
vided to the FC is built by accumulating randomly weighted 
sensors readings along a network path. Energy efficiency is 
pursued if the number of nodes contributing to each projection 
is low. 

In [22] a peculiar form of the CS sensing matrix is proved 
to exhibit good reconstruction properties while still being 
able to reduce the number of inter-sensors transmissions. The 
structure of the sensing matrix, originaUy designed for WSNs 
with chain topology is viable of an extension to more complex 
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scenarios, provided suitably tree-structured routing paths are 
designed from the exterior of the network towards the FC. 

Recently, in 1 2 1 1 a RS based sensor network framework for 
underwater systems has been introduced. Differently from the 
works in [17], [18], [20], a single hop network is considered, 
where the sensors directly communicate with the FC. Every 
sampling interval 6t, each sensor senses the field and transmits 
it directly to the FC with probabiUty p. The sensing probability 
p is suitably chosen in order to let the FC acquire sufficiently 
many data for field reconstruction. The approach in [21] favor- 
ably merges the RS procedure with a random access protocol, 
thus obtaining a significant reduction in both the consumed 
energy and the occupied bandwidth. 

The main reasons why the above discussed methods achieve 
significant performances in efficient use of network resources 
is either in the fact that, during a sampling interval, only a 
sub-set of the sensors measures the field by means of RS, or 
in the fact that CS acquisition is actually realized jointly with 
the routing procedure. 

Recent studies [20], [22], suggest that in case of spatially 
sparse signals the energy/bandwidth of RS and CS efficiency 
deteriorates, since reconstruction accuracy is guaranteed only 
when a large number of sensors contribute to the measure- 
ments. In [20], the authors remark how it is difficult to design a 
RS matrix suited to sparse signals and still allowing an efficient 
network routing. In [22] it is pointed out that, on spatially 
sparse signals, CS techniques still guarantee reconstruction 
accuracy but at an increased number of measurements (e.g. 
M up to 50% of N) whereas, for the same values of M and 
N, RS only opportunistically achieves reconstruction. 

In several WSN applications, the sensed field indeed con- 
tains local fluctuations and abnormal readings, and it is well 
modeled as a spatially sparse signal. This motivated us to 
concern ourselves with the design of a sensing matrix suited to 
spatially sparse signals, as described in the following Section. 
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Fig. 2. Bidimensional sequence z\n\,n2\: column-wise accumulation for 
computation of the measurements j/f") [m] 



V. CS USING Radon-like random projections 

In a nutshell, we aim at devising a sensing matrix that repre- 
sents a spatially sparse field z in a domain such that dropping 



N — M components does not prevent signal reconstruction. 

Recalhng that the Radon Transform has the dual properties 
of i) compressing spatial domain straight lines into transform 
domain pulses [23], ii) expanding spatial discrete pulses into as 
many nonzero values as the number of considered Radon pro- 
jections [24], we recognize that the Radon Transform provides 
a mean for redundant representation of sparse signals built 
by spatially isolated pulses. Thereby, we resort to a spatially 
sparse signal CS scheme inspired to the Radon projections 
computation. 

To elaborate, let us consider the bidimensional sequence 
z[ni,n2] of finite size A^i x A7'2and let us present few 
examples of measurements computed in analogy to the Radon 
projections. 

First, let us consider the colunm-wise accumulation of a 
randomly weighted version of z[ni,n2\: 

N2-1 

y(o)[m] = ^ [i,m], m = 0,...,N,-l (6) 

i=0 

where <^m^[«],« = 0, . . . A^i — 1, m = 0,...,N2 — 1 are 
i.i.d. zero mean Gaussian random variables of equal variance 
cr^, i.e. ipm\i] ~ AfiO.a'^); an example of the formation 
of the measurements y(")[m] is illustrated in Fig. 2. The 
measurements y^^^ [m] differ from a Radon projection in that 
each sample is randomly weighted. Besides, they differ from 
classical CS measurements in that each measure is obtained 
only from a subset (namely a column) of the values in 
z [ni,n2] rather than on all the samples z [711,712]. 

Definition of the row-wise random projections of z [rii , 712] 
is straightforward, namely: 



ATl-l 

E 

i=0 



(7r/2) 



\i]-z[m,i], m = 0,...,N2-l (7) 



By analogy, we can define the diagonal-wise projection of 

z [711,712] as well: 

( Ni-l-m 

^ ^(Z^^^[i]z[i + m,i], 

i=0 

m=0, A^i-l 

JV2-l-|m| 

^ ^<f;:/%]z[i,i + m], 



.,(V4)r-ml = 



[m]={ 



(8) 



i=0 



771 = 



-No 



1, 



-1 



The expressions (6)-(8) represent randomly weighted accu- 
mulations of z [rii, 712] over ridge paths. Let us now generalize 
the above expressions by regarding them as obtained by 
column-wise accumulation of a suitably rotated version of the 
sequence 2 [711, 772]- 

Let z^'^"^ [711, 712] be a 7?p-radiant clockwise rotated^ version 
of the image z [ni, 772]. The size of 2;(''r') [ni, 772] varies with 

^Formally, we obtain the iJp-radiant clockwise rotated version of the image 
z [ni, 712] by regular sampling of the rotated field 

s^'^"^ {x,y; t) = s^^"^ (a;cosi?p + 2/ sini?p , a; sini?p — ycosiJp; t) . 
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■!?p, and we denote as K{p),p = 0, ... P — 1 the number of 
columns of z^^"^ [ni,n2]. 

We generalize the definitions (6)-(8) by considering the 
collection of the random projections of z^'^''^ [Tii,n2\ over a 
finite set of P directions 



(9) 



for m = 0, . . . K{p) — 1, p ~ 0, . . . P — 1. Since the accumu- 
lation in (9) recalls the column-wise accumulations employed 
in the computation of the discrete Radon transform, we refer 
to the projections in (9) as Radon-like random projections. 
Let us now collect the P Radon-like random projections: 

y(^p) = [j^(^p) [0] . . . yi^p) [K{p)-l]Y,p = 0,...P-l 

in a measurement vector y: 

y=[y(''o)^..y(^P-.)Y 

The measurements y are computed from z as: 

y = (10) 

using the Ni x A'2 random sampling matrix defined as: 



('?o) 



where = 0, ... P — 1 are the suitably defined sparse 

random matrix realizing the accumulation in (7)-(8). For the 
sake of clarity, we report the matrices and 
corresponding to the horizontal and the vertical projection, in 
(5) and (4) respectively. 

Let us remark that, even though the samples contributing to 
each y component are the same that would have contribute to 
a specific Radon projection of z, due to the random weighting 
the measurements vector y is not -and it is not even easily 
related to- the Radon transform of z. 

In the foUowings, we demonstrate that the conditions given 
by the CS theory for reconstructing z from y stand. Before 
turning to mathematics, let us observe that if the image 
z [n\ , 7x2] is built by sparse isolated pulses, each pattern 
contributes, a part a suitable randomly weighting, to each 
one of the P Radon-like random projections y'^^''^[m], m = 
0, . . . K{p) -1, p = 0,...P-l. Thereby, the set of P Radon- 
hke random projections is a P-redundant representation of 



the pulse. This motivated us to formally demonstrate that the 
Radon-like random projections y'^'^''^[m\, m = 0, . . . K{p) — 
1, p = 0,...P— 1 satisfy the conditions of a CS measurement 
set. 

A. Restricted Isometry Property of the Radon-like sampling 
matrix 

To formally state that the above introduced simplified 
samphng structure is feasible for accurate field reconstruc- 
tion, we shall demonstrate that the random measurements y 
evaluated as n (10) are consistent CS measurements, i.e. they 
substantially preserve the information of the sampled sequence 
z[n-i,n2\. 

Formally, we need to prove that the sampling matrix 
satisfies the condition known as Restricted Isometry Property 
(RIP). Specifically, let us denote by £2 =^||z|p the quadratic 
norm of the vector z. If the matrix satisfies: 

il-5)£,<\\^Rz\f <il + 5)£, (11) 

with high probability, then, any sparse sequence z can be 
perfectly recovered from CS measurements y = $i{Z [13]. 
The RIP in (11) asserts that the measurements energy 
=^||$jiz|p is strongly concentrated around the value Sz- 
Prehminary results on the RIP property of a Radon-like 
CS matrix appear in [25 J. In Appendix A, we extend these 
results and, following the approach in [30], we show that the 
following property stands. 

Property 1: If the nonzero entries in are i.i.d. zero mean 
Gaussian random variables with equal variance cr^ = 1 /P, the 
following concentration inequality stands: 

Ml^y - £z\ > S} < e (12) 

provided that P > 2K'^Cilog{2/e)/S^ being C2 a suitable 
constant. From Property 1, the RIP property in (11) immedi- 
ately follows. 

B. Further Remarks 

The samphng matrix i>ij differs from the full samphng 
matrices usually adopted in CS since it is sparse. In the 
following, we show how the sparsity of the Radon-like sensing 
matrix can be leveraged on to significantly simpUfy the 
CS measurements computation in a WSN, so as to reduce the 
consumed energy and employed bandwidth. 

The definition of y'^^^\m] in (9) given above is con- 
sistent and useful from an analytical point of view. When 
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turning to evaluation of y*^''p^[m] in a WSN, the evaluation 
of z'^'^"^ [ni,n2\ is not accomplished, and the measurements 
computation is realized throughout the data gathering stage. 

Finally, the concentration inequality (12) guarantees that 
the sequence z can be recovered from the measurements y 
with high probability, as far as the number P of considered 
projections increases. Since CS convergence is assured only 
in probabihty, the CS measurement experiment could be 
repeated. In a WSN, data are periodically sensed and routed 
to the FC, and a small probability of mis-reconstruction can 
be tolerated since it can be recovered in the subsequent 
sampling interval. Furthermore, integration of independently 
drawn measurements acquired in a WSN during different 
temporal intervals can be performed at the FC. This interesting 
research issue is deferred to further studies. 

VI. Radon-like CS in a WSN 

In a WSN application scenario, the P Radon-hke projec- 
tions correspond to randomly weighted sums of the values 
sensed by different sub-sets of WSN sensors. The sums can 
be evaluated using different techniques. 

Following the approach in [19], the sensors within a subset 
can synchronously transmit their weighted sensed values, and 
the sum can be reaUzed at the FC by on air analogical 
superimposition of received signals. This protocol requires 
a strict control of the power received by the FC from each 
sensor. Precisely, each sensor node needs to estimate the 
channel seen towards the FC, in order to pre-compensate 
the transmitted value according to the channel attenuation. 
Thereby, although feasible in principle, this approach requires 
a sophisticated processing and a tight power control by the 
sensor nodes. 

According to a data gathering paradigm, the projections are 
computed within the network by sub-set of sensors while they 
are forwarding their sensed values to the FC. The sums in (9) 
can then be reaUzed by routing and accumulating values of 
z [ni , 712] over suitably tilted paths in the network grid discrete 
support. Here we refer to such a data gathering approach, and 
we infer some consequences from the pecuUar sparse structure 
of the matrix on the computation procedure. 

Firstly, we observe that in every row, the non-zero coeffi- 
cients of the matrix are arranged so as to obtain y^'^p^ [m] 
as the sum of the values of a column of the rotated image 
2;(''p). When collecting the measurement in a WSN, each 
projection y^'^''^[m] can then be computed by accumulating 
measurements throughout a specific, suitably tilted, grid path. 
Secondly, we observe that in every column of the matrix 
there are only P non-zero values. Hence, each value z [i,j] 
shall contribute only to P out of M projections y'^'^P'' [m]. 
When realizing the Radon-like CS in the sensor network grid, 
each sensor shall transmit its value P times. 

Based on these premises, we recognize that the sparse 
structure of the matrix results into two main features of 
Radon-Uke projections computation in a WSN: 

• the computation of each projection j/'^p^ [m] is performed 
in a distributed way within the WSN, and it requires 
signaling among grid sensors which are adjacent along 
a WSN path; 



• the accumulation along different paths can then be re- 
aUzed in parallel, provided that the distance between 
contemporaneously signaling nodes is kept large enough. 

We now turn to quantifying the energy consumption and 
bandwidth occupancy required for realizing the Radon-like CS 
by means of a data gathering procedure in a WSN. In VI- A 
we evaluate the energy consumption and bandwidth occupancy 
of the Radon-like CS scheme in a WSN. Next, we compare 
these results with selected state of the art schemes, namely 
those referring to the conventional network and to the Random 
Sensing approaches described in Sect.I, respectively detailed 
in subsections VI-B and VI-C. 



A. Radon-Like CS efficiency 

We now evaluate the allocated bandwidth and consumed 
energy for entirely collecting the measurements in a time Tc 
in the WSN scenario in Sec.II. The actual bandwidth occu- 
pation and energy consumption depend not only on the WSN 
structure but also on the adopted data gathering procedure. 
Without loss of generahty, we refer to the suboptimal data 
gathering algorithm sketched out in Appendix B, and here 
briefly summarized for the reader's convenience. 

According to the algorithm in Appendix B, the sensors 
transmit their readings to the FC by data gathering through 
suitable multi-hop paths. Figs. 3-4 illustrate, as an example, the 
multi-hop paths selected for the computation of the projections 
y('r/2)jj^j^y(7r/4)j^j yyifjuji ^ nctwork quadrant. 
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Fig. 3. Data gathering algorithm: spatial organization of the sensors' 
transmissions for horizontal projections evaluation in the first quadrant. Time 
stamps indicate when the node starts transmitting. 

A deterministic (collision free) Time Division Multiple 
Access (TDMA) is adopted. Signaling takes place between 
adjacent nodes only. Parallel transmission of sufficientiy apart 
nodes are considered. 

In Appendix B, we evaluate two parameters that directly 
affect the energy and bandwidth consumption of the data 
gathering algorithm, namely total number of single hop trans- 
missions A'tx. and the total number A'xs of time slots needed 
to collect all the sensors' readings to the FC. 
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Overall, we can express the total consumed energy^ for the 
Radon-Like CS based approach as 



Fig. 4. Data gathering algorithm: spatial organization of the sensors' 
transmissions for diagonal projections evaluation in the first quadrant. Solid 
arrows indicate pi) paths, dashed arrows indicate p2) paths. Time stamps 
indicate when the node starts transmitting. 



(15) 



being the overall number of transmission. 

Besides, in order to assure that the measurements are 
collected within a maximum time of Tc, the packet time Tp^^ 
needs to meet 



^(RL) 
P 



T 

c 



so that we have 



d{P} 



(16) 



The occupied bandwidth, defined as the packet length in 
bits L over the packet transmission time, for the Radon-Like 
approach is finally evaluated as 



(17) 



The total number of node transmission A'tx for the al- 
gorithm under concem comes out to vary linearly with the 
network size: 

iVxx « 7p . AT, (13) 

being 7p a scalar factor depending on the adopted P pro- 
jections. The guidelines for calculating 7p are given in the 
Appendix B where we also evaluate the value of 7p for 
different sets of directions dp. 

Based on the same data gathering algorithm, we have 
evaluated the number of time slots needed to collect all the 
sensors' readings to the FC. The number of time slots tums 
out to vary just with the square root of the network size, that 
is: 

Nts « Sp ■ Vn, (14) 

being 5p a factor depending on the considered projection 
directions 'dp. Appendix B reports the guideUnes for the 
evaluation of the parameter 5p for different sets of directions 

dp- 

We observe that, as a consequence of the sparsity of the 
Radon-Uke sensing matrix, i) the number of transmission 
varies only linearly with the network size N\ n) being the 
algorithm paralleUzed, the number of time slots varies linearly 
with the square root of the network size N. 

With these results, we are able to evaluate the allocated 
bandwidth and consumed energy for entirely collecting the 
measurements in a time Tc in the WSN scenario described in 
Sect.n. 

Projection evaluation according to the data gathering 
scheme detailed in Appendix B accounts for a series of 
transmissions among neighboring nodes. The energy spent 
for a single-hop transmission is given by i?sH = Gd^y^T^^^ 
being dsn = d ox \f2d the distance for horizontally, vertically 
and diagonally adjacent nodes, and the time needed for 
packet transmission. 



B. Conventional network efficiency 

As a benchmark for reconstruction quality and en- 
ergy/bandwidth consumption comparison purposes, we refer 
to a conventional network in which, at time t, all of the sensors 
sample the underlying field and transmit their measurements 
via a TDMA scheme. iV time slots are organized in a frame 
structure of duration = NTp. Recalling that the FC needs 
to reconstruct the sensed field every it must be Tf < Tc- 

In the conventional network, all the sensors access the 
medium for packet transmission in a TDMA fashion. Thus, 
the conventional network packet transmission time Tp'"'^^'' , is 
given by T^^"^^ = Tc/N. The overall energy consumed by 
the conventional network is then evaluated as 

^coN = E P^'^rf ON) = ^ G 4Tf ON) 

k k 

Since the FC is located at the center of a regular grid, 
as depicted in Fig.l, the value of dk can be expressed in 
terms of the inter-sensor distance d: dk — d + 
where i, j range in {-[Ni/2\.[Ni/2\),{-[N2/2\, L^2/2J), 
respectively; then, ii^coN rewrites as 

i j 

The sums over i and j in (18) can be approximated as 

i j ^ ' 

Denoting Ni '^=aiy/N, N2 = a2VN, the overall energy con- 
sumption for the conventional network sums up to 

E^o^ = ^^^^(^^ + ^1) G d^iV^Tf ON) (19) 

'We herein approximate also the inter-sensor distance in the diagonal 
direction to d. 
''Recall that J2?=i = + l)(2n + l)/6. 
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By setting T^'^°^^ = T^/N in (19) we come up with: 



CON 



aia2(ai + ai) 

48 



G d^NTr 



(20) 



The occupied bandwidth for the conventional network is: 

(21) 



L L 
-BcoN = — = 



■,(CON) 



C. RS Efficiency 

We now summarize the energy and bandwidth consumption 
of the approaches proposed in [21]. Therein, a RS procedure, 
allowing only a randomly chosen sub-set of sensor to acquire 
the measurement, is coupled with both a TDMA scheme and 
a random access scheme. 

In the RS/deterministic access (RD) scheme, the FC ran- 
domly chooses a set of M sensors, M being a sufficient 
number of measurements for satisfactory field reconstructions 
and it broadcasts the addresses of eliged nodes through the 
network. The selected nodes acquire the measurements and 
transmit their readings to the FC via a TDMA deterministic 
access scheme. 

As in this scheme only M nodes need to share the TDMA 
frame, the packet transmission time is Tp^^^ ~ Tc/M. 
Consequently, the occupied bandwidth for the RD scheme is 



Brd = M— 



(22) 



In order to evaluate the consumed energy, let us consider the 
set Cm;N collecting all the possible configurations of M out 
of N nodes. The energy consumption of a given configuration 
c € Cm;N of M is expressed as 



where the sum over the index kc spans the M sensors within 
the combination c. Thereby, the energy consumption of each 
combination depends on the distances of the M nodes from 
the FC. In this respective, we evaluate the energy consumption 
of the RD scheme as the average over all the possible 
combinations: 



Erd — 



1 



K 



M;N 



(23) 



where Km-n = (^) cardinality of Cm;N- After 

rewriting the above equation as 

Erd 



Km:N 



we recognize that in the overall sum over the Km-.n combi- 
nations, the energy spent by each and every network sensor 
appears in (^Z^) terms, corresponding to the combinations it 
belongs to. Therefore, the above sum can be rewritten as 



Erd — 



K 



M;N 



N -I 
M-1 



fe=i 



(24) 



Following the same steps already detailed for the conventional 
network we have 



N 48 
which, given that Tp^^^ = T^M, reads as: 

■"^ " 48 ■= 



(25) 



(26) 



The energy consumption and occupied bandwidth perfor- 
mance of the scheme in [21] need to be addressed in a sUghtly 
different way w.r.t. the previous cases. CS theoretic results 
determine the number M of measurements needed at the FC 
to correctly restore the sensed field; within an observation 
time, this constraint corresponds to a required percentage qs 
of correctly received samples at the FC. Because of possible 
collisions, the required percentage Qs needed at the FC does 
not translate straightforwardly into a sensing probability ps- 

In [21], the authors estabhsh the following relationship 
among the sensing probability Ps and the probability qs of 
correct packet reception: 



: Ps exp -p, 



2LN 



(27) 



being B the available bandwidth. Remarkably, given qs, a 
minimum bandwidth is required in order to assure that a 
feasible value of Ps exists. Thereby, if the available bandwidth 
is not accurately dimensioned, small values of ps do not 
provide enough measurements at the FC, whereas large values 
of Ps cause too many collisions. 

In [21], the authors provide an expression, standing for 
small values of qs, of the minimum bandwidth as a function 
of the desired qs. Following the guidelines in [21], we have 
extended such result to accommodate for large values of qs, 
too. Specifically, we came up with the relation: 



{2Ne-qs + l) if g« < e"! 



L f 2N 



To \-lnqs 



(28) 



+ 1 if >e-i 



The RS settings are therefore assigned as follows. Firstly, 
the desired qs is fixed according to the reconstruction quality 
constraints; secondly, the minimum needed bandwidth is eval- 
uated according to (28); finally, the selected bandwidth value 
is employed to derive the needed sensing probability Ps using 
(27). 

With these positions, the packet transmission time Tp is 
determined from the employed bandwidth as t},^^^L/B^^. 
Besides, the energy consumption is determined by the value 
of Ps', the average network consumed energy is evaluated as 

i^RR = E|^p,pWTf'^)| = |Ar2T('^) (29) 
By expUcating the value of Tp^\ the expression in (29), 
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depending on the value of qs, reads as follows: 



Err = < 



T,- 

24 {2Neqs + l) 
PsN^ Tc 



24 



2N 



+ 1 



if qs<e ^ 
if qs > 



(30) 



For high values of N, the expressions in (30), can be globally 
approximated as 



Err ~ —NTc 
48 



(31) 



D. Further remarks 

Before turning to the numerical performance evaluation, a 
few remarks are in order 

The above analysis has pointed out that the consumed 
energy and allocated bandwidth adopting the Radon-Uke CS 
scheme grow only with the square root \/N of the network 
size, where as those of selected state-of-the-art approaches 
vary with the power N. Therefore, the Radon-like CS scheme 
yields a reduced energy consumption for each node, as well 
as a parsimonious bandwidth use for collecting data over all 
the grid. The gain is more and more evident as the network 
size (i.e. the covered area) increases. 

The actual gain in terms of energy and bandwidth depends 
on the constants 7p,5p which grow with the number of 
considered projections. The advantages of the Radon-like CS 
scheme are expected to be evident on spatially sparse signals, 
where low values of P (e.g. P = 3) enable reconstruction, 
whereas RS data gathering algorithm require a high percentage 
of the samples to reach the FC in order to achieve satisfying 
reconstruction results [22]. 

As far as the medium access scheme are concemed, the 
herein devised Radon-hke CS scheme is reahzed via a TDMA 
technique, just as the RD scheme; therefore, it implies an effort 
of synchronization and scheduling. Nonetheless, different data 
gathering procedures can be envisaged realizing the Radon- 
Uke CS using a random access criterion; this issue is left for 
further studies. 

Finally, these results depend on the pecuhar structure of 
the Radon-like sensing matrix and although derived for 
a particular data gathering algorithm, can be generalized to 
different Radon-hke CS measurements computation schemes. 



Vll. Numerical simulations 

In this Section, we analyze the performance of our Radon- 
like CS scheme both in terms of reconstruction accuracy and 
of employed network resources. Firstly, in Subsect.Vll-A we 
present a few results showing that an accurate reconstruction 
of a spatially sparse signal can be achieved by the RCSM 
using a feasible number of Radon-like projections. Secondly, 
in Subsect.VII-B we demonstrate the energy and bandwidth 
gain achievable when adopting the Radon-Uke CS scheme in 
a WSN. 



A. Radon-like CS reconstruction accuracy 

With reference to the network model in Sect.Il, we consider 
a WSN made up by a square grid of = 64 x 64 = 4096 
sensors. The sensed field z[ni, ^2] is built up by 7 repetitions 
of an elementary 5x5 pattern, differently scaled by factors in 
the range (0.5 — 0.85]; Fig.5-a) describes an example of the 
field z[ni,n2]. This field adhere to the so called pulse stream 
signal model described in [29], and the therein presented 
Algorithm 2 is employed for computing the reconstructed field 
z[ni,n2]. 

We have first evaluated the reconstruction accuracy of 
the Radon Like CS scheme, using different number P of 
projections, under the assumption of noise free observations. 
Specifically, we have tested the reconstruction accuracy when 
only projections along rows and columns of the network grid 
are considered (P = 2, M = 128), when projections along the 
rows, the columns and along the tt/A oriented diagonal of the 
grid are considered (P = 3, M = 255) and finally when also 
projections along the 37r/4 oriented diagonal are considered 
(P = 4,M =^ 382). 

Tab.I reports the MSB of the reconstructed field z, i.e. 
MSE = (z — z)'^(z — z), achieved after 20 iterations of the 
reconstruction Algorithm 2 and averaged over 10 runs corre- 
sponding to different pulse locations. Results in Tab.I show 
the effectiveness of the Radon Like CS scheme in sensing a 
sparse field with a restrained number of measurements. 



M 


128 


255 


382 


P 


2 


3 


4 


■&V 


0,7r/2 


i?p = 0,7r/2,7r/4 


i?p = 0,7r/2,±7r/4 


Radon Like CS 


0.0056 


0.0019 


0.0011 


RS 


13.3 


7.604 


4.721 



TABLE I 

Reconstruction ACCURACY (MSE) obtained by the Radon-Like 
CS scheme and by the RS scheme (Af = 64 X 64 = 4096). 



To visually assess the Radon Like CS scheme reconstruction 
accuracy, we show in Fig.5-b) the reconstructed field z[n\, 722] 
obtained by measuring the field in Fig.5-a) with P = 3 
projections. 

For the sake of comparison we have also evaluated the 
reconstruction accuracy obtained by the RS [21] scheme under 
the same experimental settings. In Tab.I we recognize that, 
for the selected range of measurements (M < 382) the RS 
doesn't allow to capture the sparse nature of the sensed field. 
To obtain the same reconstruction accuracy of the Radon Like 
CS scheme, namely a MSE equal to 0.0036, the RS requires 
be run with a number of measurements M ~ 2500 « 50%A'^. 

Similar results have been obtained by scaUng the number 
of network nodes to A'^ = 80 x 80 = 6400 and considering a 
sensed field z[ni. 71,2] composed by 8 repetitions of the elemen- 
tary pulse. The MSE of the Radon-like CS scheme, averaged 
over 10 runs, is reported in Tab.II. Also in this case the MSE 
obtained by the Radon Like CS scheme with P = {2,3.4}, 
corresponding to M = {160, 319, 478} measurements, proves 
to exhibit satisfactory reconstruction quaUty. For comparison 
sake, we observe that in these experiments, M = 3500 were 
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a) b) 

Fig. 5. Original field 2[ni,n2] (a) and field i[ni,n2](b) after sensing 
performed according to the Radon Like CS with P = 3 projections. 



needed by the RS to achieve analogous performance, namely 
a MSE equal to 0.006231. 



M 


160 


319 


478 


P 


2 


3 


4 




0,7r/2 


i9p = 0,7r/2,7r/4 


■dj, = 0,7r/2, ±7r/4 


Radon Like CS 


0.0045 


0.00152 


0.000901 


RS 


8.602 


6.472 


4.506 



TABLE II 

Reconstruction accuracy (MSE) obtained by the Radon Like 

CS SCHEME AND BY THE RS SCHEME (AT = 80 X 80 = 6400). 



Finally, we have tested the Radon Like reconstruction accu- 
racy when noisy acquisition are considered so that the CS 
measurements can be modeled as: 

y = $z + n 

where the vector n gathers samples of white zero mean 
Gaussian noise with variance cr^j. In Tab.III we report the 
reconstruction accuracy obtained when acquiring with N = 
64 X 64 — 4096 sensors a field composed by 7 patterns 
via the Radon Like CS scheme for different values of cr^. 
Results in Tab.III show how the presence of noise in the 



M 


128 


255 


382 


P 


2 


3 


4 




0,7r/2 


■dp = 0,7r/2,7r/4 


i?p = 0,7r/2,±7r/4 


al = 0.5 


0.0065 


0.0034 


0.0033 


< = 0.7 


0.007 


0.0039 


0.0036 



TABLE III 

Reconstruction accuracy (MSE) obtained by the Radon Like 
CS; noisy measurements (iV = 80 X 80 = 6400). 



acquisition process does not severely affect the reconstruction 
performance of the Radon Like CS approach. 

To recap, the above results show that, as a rule of thumb. 
Radon Like CS requires M sa P\/N measurements for 
representing a spatially sparse field, whereas the RS requires 
a large percentage of the measurements A/_rs ~ aN to be 
correctly received {e.g., P = 3 and a = 50% in the above 
experiments). Overall, the Radon Like CS scheme allows 
sensing and reconstructing of a spatially sparse field with far 
less measurement w.rt. to state of the art techniques such as 
the RS presented in [21]. 



B. Radon Like CS efficiency 

We now show that, besides using a restrained number of 
measurements, the Radon Like CS presents significant advan- 
tages in terms of energy and bandwidth needed to disseminate 
sensors readings towards the FC. 

We evaluate the performance of Radon Like CS data gath- 
ering scheme (RL) by evaluating the consumed energy and 
occupied bandwidth under the following assumptions 
al) the number of measurements is chosen large enough to 

yield a satisfactory reconstruction accuracy, quantified by 

a MSE < 10-3; 
a2) the coherence time of the sensed field is fixed to Tc = 

2500s; 

a3) the dimension of the transmitted packet is set to L = 
1Kb. 

With reference to the experimental setting described in 
Sect.VII-A, condition al) implies P = 3 projections. 

We start by computing the occupied bandwidth required by 
the RL scheme to guarantee correct sensing and reconstruction 
of a spatially sparse field. Fig. 6 plots the occupied bandwidth 
evaluated according to (17) vs the network size in terms of 
number of sensors. 

For the sake of comparison, in the same figure we also 
report the performance of the RS, implemented both using 
a deterministic access scheme (RD in the legend) and the 
random access scheme (RR in the legend) described in [21] 
and analyzed in Sect.VI-B; besides, we report the results for 
the conventional network (DD in the legend) described in 
Sect.VI-B. The occupied bandwidth is computed according 
to (21), (22) and (28) for the RL, DD, RD and RR scheme 
respectively, under the same assumptions al)-a3). These set- 
tings imply ps = 0.5 for the RD scheme and Qs = 0.5 
for the RR one, whereas the conventional network, bringing 
all the measurements to the FC, always achieves perfect 
reconstruction. 

Occupied Bandwidth 




1000 1500 2000 2500 3000 3500 4000 



Network Nodes 



Fig. 6. Bandwidth occupancy versus number of nodes N. 

The bandwidth employed for the Radon-like data gathering 
scheme is significantly reduced, not only w.rt. the conven- 
tional network, but also w.rt. the RS, both using deterministic 
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Occupied Bandwidth 
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Fig. 7. Bandwidth occupancy versus number of nodes N. 
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Fig. 8. Energy consumption versus number of nodes N. 



Bandwidth/Energy Pairs 
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Fig. 9. Energy/Bandwidth pairs. 



and random access. The above results, obtained by considering 
different number of measurements and equal reconstruction 
accuracy for the three algorithms, are explained by the effi- 
ciency of the Radon-like data collection algorithm, exploiting 
only single hop, parallel data transmission. 

We have then pushed further the comparison to find under 
which conditions the Radon-like and the RS schemes equally 
perform. In doing so, we have evaluated the bandwidth oc- 
cupied by the RD and RR schemes when only 10% of the 
sensor readings are required to be correctly received by the 
FC. Fig. 7 reports these results showing that the RL scheme 
still favorably compares with the RR scheme, while requiring 
just slightly higher bandwidth w.r.t. the RD scheme. Therefore, 
the RL scheme overcomes the RD and RR ones unless they 
use a very low percentage of the sensors' readings for field 
reconstruction. 

Let us now compare the energy consumption under the 
assumption that the occupied bandwidth, and consequently the 
packet transmission time Tp = LjB, are the same for the 
different schemes. 

Let us remark that each scheme requires a different number 
of transmissions in order to let the FC acquire the needed 
measurements. Thereby, the overall process of sensing is 
accomplished in different time intervals by the different data 
gathering techniques. In these experiments, we have set Tp = 
0.61 so that also for the DD scheme, comprising the maximum 
number of packet transmission, the data gathering time meets 
the constraint on the field coherence time (namely ISSTp < Tc) 
even for = 4000. Fig. 8 reports the energy consumed for 
different number of network sensors N; we recognize that the 
employment of the RL scheme drastically reduces the energy 
consumed by the data gathering algorithm. 

Finally, we are interested in comparing the resource em- 
ployed by the different schemes to accomplish the sensing 
process exactly in the same time. To reach this goal, each 
scheme occupies a different occupied bandwidth. In Fig. 9 we 
illustrate the bandwidth-energy scatterplot of the RL, RR, RD 
and DD schemes for different network sizes. For the RD and 
RR we have considered two cases, that is when 50% and 30% 
of the sensors readings are required at the FC for satisfactory 
field reconstruction. 

The energy/bandwidth pairs draw different trajectories while 
the number of nodes N increases. We clearly recognize a 
systematic energy and bandwidth saving of the RL scheme 
w.r.t. to the competitors. Moreover, both the consumed energy 
and the occupied bandwidth exhibit far smaller variations with 
the number of network nodes than what happens with the RD, 
RR and DD schemes, making the RL approach fully scalable 
in terms of network nodes. 

As a final remark, we observe that energy and bandwidth 
gain yielded by the Radon-like approach is directly based 
upon the favorable matching between the Radon-like matrix 
structure and the spatially sparse structure of the sensed field. 
Based on these results, we envisage a twofold extension of 
our work, that is 

• in accounting an irregular sampling grid; 

• in exploiting non straight sampling path. 
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These conditions may be encountered, for instance, in ve- 
hicular networks or citizen sensing networks [33], where the 
disposition of the nodes are far from being regular, and the 
sampling path should adapt to the routing paths, which in turn 
basically depend on the streets and buildings disposition. 

To sum up, we consider this contribution as an intermediate 
step towards the finding of a general relation between the 
compressive sensing of a finite innovation rate signal and its 
realization by efficient routing algorithms in a realistic WSN 
scenario. 

VIII. Conclusion 

In this paper, we have addressed the efficient compressive 
sampling of spatially sparse signals in sensor networks. Specif- 
ically, we have introduced a peculiar CS samphng scheme 
for spatially sparse bidimensional signals. We have analyti- 
cally demonstrated that our scheme satisfies the theoretical 
conditions required for CS signal reconstruction. Then, after 
devising a distributed data gathering scheme for collecting 
the CS measurements in a WSN, we have characterized the 
scheme both in terms of consumed transmission energy and 
occupied bandwidth. The scheme outperforms state-of-the- 
art schemes for spatially sparse fields, and it represents an 
intermediate step towards the definition of routing procedures 
well-suited to the characteristics of the signal a realistic sensor 
network is faced to. 

Appendix A. RIP property of the Radon-like 

MEASUREMENTS MATRIX. 

Property: If the nonzero entries in are i.i.d. zero mean 

Gaussian random variables with equal variance cr^ = l/P, 
the following concentration inequality stands: 

^T{\Ey-£,\>5}<e (I.l) 

provided that P > 2K^Cilog{2/e)/d^ being C2 a suitable 
constant. 

Demonstration: 

The condition in (I.l) can be demonstrated as follows. 
Let us consider the sample energy £y of the measurements. 
Being £y a sample moment, we here invoke its asymptotical 
normal distributions [31]. Although this hypothesis is not 
necessary for the RIP to stand, it allows us to straightforwardly 
evaluate the minimal number of projections P required for 
if-sparse field reconstruction, and therefore we retain it in the 
foUowings. 

The Chernov bound [32] for a normal random variable 
establishes that the probability that the random variable differs 
from its mean is limited by a term exponentially decaying 
with its variance. By application of the Chemov bound to the 
random variate £„, we obtain: 



PT{\£y-E{£y}\ >S}< 2-cxp 



Thereby, in order to demonstrate (12) it suffices to show that 
E {£y} = £z and Var {£y] = o{P) for the case under concem. 
By definition we have 



p-l K{p)-1 
p=0 m=0 



Let us denote by 

= E W^*''''' [i, m],p = . . . p-l, m = 0, . . . K{p)-1 

i 

We recognize that are i.i.d. normal random variables, 

with zero mean and variance equal to 

i 

With these positions, we recognize that 

P-l K(p)-1 

E{f.}=E E Vary^) 

33=0 m=0 

P-l 



p=0 i 



As far as the variance of £y is concerned, we obtain 

p-l K(p)-1 



V^r£y=Y: E Var(yCf.)) 

33=0 m=0 

P-lK(p)-l / , , , , 2\ 

= 1: E (i^ {(£'■) }) 

Since the variables ym are zero-mean normally distributed, 
their fourth-order moments satisfy 

E{(yL^^))'}=3-(VareO' 



so that we obtain 



p-l A'(p)-1 



Var£:, = 2^ (Varj/^f.)) 

p—O rn—Q 

Observing that, for a K-sparse signal, the maximum value of 
{^^^"^['■2^^]) is achieved in case of K aligned pulses, 
we recognize that the following inequaUty stands: 

P-iK(p)-i 

(Vary(;f^)) =<E E El-^'^^lH,-]) 

i-2 

P-\K(ji)-\ 

^<E E El-^'^^i^i'-]) -^c- 

p=0 m— i\ 

= a%P ■ • KC2 

being C2 = max(„^_„2) ^ ["i) '^2]^- 

Finally, we obtain that the variance Yax{£y} is upper 
bounded by: 

(1.2) 



1 



\'^r{£y}<--K^Cl 

In (1.2) we recognize that the variance Varjfj,} decays 
as 1/P. Furthermore, by comparison of (1.2) and (12) we 
recognize that the RIP is verified provided that the number 
of projections P satisfies: 

P> 2i^2c|log(2/e)/52 
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Appendix B. Within WSN Radon-like projections' 

COMPUTATION 

Let us consider a regular network composed by iV = N1N2 
sensors as in Fig.l where, without loss of generality, we 
assume that iVi and N2 are odd valued, i.e. Ni = 2Ni + 1, 
N2 = 2N2 + 1, so as to identify a central column where the 
FC is located. 

We discuss a simple sub-optimal procedure to collect all the 
projections to the FC using a TDMA access scheme. 

Let us subdivide the network in 4 quadrants, and let us 
consider first the horizontal projections y^^/^'>[m], obtained 
by accumulating randomly weighted values along the network 
rows. In each quadrant, the data gathering process starts at 
the outer nodes. The external node in each row measures the 
field, computes the product of the reading with a randomly 
selected coefficient, encodes this value in a packet of L bits, 
and transmits it to the neighboring node in the horizontal 
direction. The neighboring node, once the packet from the 
outer node has been received, measures the field, multiplies 
the reading by the random coefficient, and sums it to the value 
received by the outer node. The overall process continues until 
the nodes in the central column are reached by the data flow 
and are then ready to propagate the projection values to the FC. 
Once the FC has received the horizontal projection values from 
the first quadrant, the same operations are serially performed 
in the remaining 3 quadrants. 

The projection values of each quadrant are therefore com- 
puted by evaluating partial sums and propagating them towards 
the nodes in the central column; then, the projection values are 
transmitted to the FC via a multihop route along the central 
column. 

Let us now evaluate the number of transmission required 
for computing the horizontal projections. For collecting the 
contributions within a network quadrant we need 

• A^i transmission to reach the central column for each one 
of the A^2 + 1 rows; 

• S/^'o ^ transmissions to propagate the projections value 
towards the FC along the central column. 

Accounting for the four quadrants, the overall number of 
transmissions for horizontal projections evaluation sums up 
to 

By denoting A^^i ^ ai^/N, N2 =^a2>//V", we can write: 

2 

N^l^^^ = (aia2 + Y)iV + (ai - a2)ViN) - 1.5 (113) 

Let us now evaluate the number of time slots in which the 
Nj'^^'^^ transmissions can be performed. Since signaling occurs 
between adjacent nodes, the propagation of the information on 
the different rows can be scheduled in parallel flows, provided 
that a suitable inter-row delay of i^q time slots is introduced 
to prevent interference among neighboring nodes. 

Let us sketch out a possible time scheduling for within quad- 
rant transmission, corresponding to the following gathering 
protocol: 



• the data gathering starts at tg = 0, on the first row of 
the quadrant, i.e. the one comprising the FC; the outest 
node transmits its randomly weighted sensed value in 
the first time slot; in the second time slot, the second 
node forwards the sum of the received data and its own, 
randomly weighted, sensed value; similarly, each node 
updates and sends the received partial sum; thereby, the 
FC retrieves the accumulation after t^^'' = NiTp, being 
Tp the duration of a time slot; 

« on the second row, the transmission begins after i^q time 
slots to avoid interference with the first row transmission; 
then, iVi time slots are needed for the partial sum to 
reach the central column and one additional time slot is 

(2) 

needed to reach the FC; the propagation ends at ^ = 
K + iVi + l)Tp; 

• on the z-th row, the transmission begins after (i — 1) • i^q 
time slots and the propagation ends at f^*"* = {i ■ i/q + 
N2 + i)Tpi 

• on the {N2 + l)-th row, transmission to the FC is 
accomplished at t'-f^^^^ = (iVst^o + Ni + N2)Tp. 

For the sake of clarity, we report in Fig. 10 a scheme sum- 
marizing the timing of nodes' transmissions when computing 
the horizontal projections [m] within a quadrant of the 

network. 




Fig. 10. Data gathering algorithm: timing of the nodes' transmissions for 
within wuadrant evaluation of y^'^^^^[m]. 

The FC is then able to collect all the horizontal projections 
in a quadrant after Ng — N21'q + Ni + N2 time slots. If the 
quadrants are visited in a serial fashion, the overall number of 
time slots to compute the horizontal projections accounts for 

N^l^^^ = i{N2i^o + Ni + N2) 

Again, by denoting A^i ai^/N, N2 =^q;2a/]V, we have: 

N^l^^^ = 2Vn [a2 + ai(l + z/q)] - (4 + 2iyo) (II.4) 

The overall protocol for evaluating [m] is illustrated 

in Fig. 3. 

An important remark is in order Despite its simplicity, this 
basic results highlights one of the major advantages of the 
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Radon-like CS scheme. Since the Radon-like projections can 
be evaluated by means of information propagation on linear 
paths in the network, the number of single-hop transmissions 
vary with the product depth and width of the network grid, 
that is with the network size N. On the other hand, since 
only single-hop transmissions are employed, the transmission 
can be parallelized and the number of time slots required for 
computation of an assigned projection set y^'^p'>[m] vary with 
the sum of depth and width of the network grid, that is with the 
square root of the network size N. This intrinsic behavior, that 
holds for different projections' directions, is the reason why 
the Radon-like CS scheme will be prove to be both energy 
and bandwidth efficient. 

With slight modifications, the above described procedure 
can be extended to the case of differently tilted paths. For 
concreteness sake, we develop in the following the calculations 
for = (vertical projections), dp = ±7r/4 (diagonal 
projections), which have been considered in the simulations 
described in this paper. 

Vertical projections y^^^ [m] can be computed using the same 
protocols adopted for the horizontal projections so that we 
have: 

AO) N2 - 1 



{2Ni + N2 + 3) 



N^l^ ^ iiNii^o + N2 + Ni) 



or, equivalently 



(0) 



{a2ai + ^)N+{a2-ai)y^{N)-1.5 (II.5) 



(11.6) 



N!f^^ =2\/iV[ai+a2(l + t^o)]-(4 + 2i/o) 



If only horizontal and vertical projections are considered 
within the Radon Like CS scheme, stemming from (II. 3) and 
(II.5) we can evaluate the overall number of transmissions 
needed to propagate data towards the FC as 



(0,7r/2) 



2aia2 



N -Z 



where we recognize how the required number of transmissions 
grows linearly with the number of network nodes 
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(0,7r/2) 



''I2N 



2aia2 



at 



For what the number of time slots needed to collect the 
sensors' readings, referring to (II. 5) and (II. 6), we have: 

^TS^"^'^ = (ai + "2) (4 + 2uo) VN-{S + 4i/o) 

where we recognize how the required number of time slots 
grows with the squared root of the network size: 

N^l'""/'^ « 62 



N 



82 = {ai + 0.2) (4 + 2i/o) 

For what the diagonal projections y^'^l^^ [m] are concerned, 
the procedure is similar to the previous case, in the sense that 
the gathering scheme is aimed at propagating the projection 
values from the outer row and column of the quadrant to the 




max{N,,N,)T 



last p2) path 



Fig. 11. Data gathering algorithm: timing of the sensors' transinissions for 
diagonal projections evaluation. 



FC by first reaching the nodes in the inner row and column. 
We serially perform the accumulations in each quadrant, as 
described for the horizontal projection, along 

pi ) all the 7r/4 oriented paths that originate from the nodes 
lying along the outer row of the quadrant and that reach the 
FC through either the central row or the central column (solid 
arrows in Fig.4); 

p2) all the tt/A oriented paths that originate from the nodes 
lying along the outer column of the quadrant and that reach 
the FC through either the central row or the central column 
(dashed arrows in Fig.4). 

The paths pi) - p2) exhibit all a length of at most 
max{Ni, N2}. Each of the A^i + A^2 + 1 projections are 
composed by at most max{7Vi, 7V2} transmissions, so that the 
overall number of needed transmissions sums up to 



N^T^ = (^1 + N2) (maxjA^i, iV2} - 1) 

= {ai + 02) max{ai, a2}N — {ai + a2)y/N 



(11.7) 



If the quadrant is processed so that we firstly perform the 
accumulations along the paths in pi) starting at to — 0, and 
secondly we perform the accumulations along the paths in p2), 
then: 

• the first projection value along the paths in pi) reaches 
the FC at ty = iVaTp; 

• the last projection value along the paths in pi) reaches 



the FC at most at t 



(Ni 



f 



= (7Vii/o + max{7Vi,7V2})Tp; 
the first projection value along the paths in 
p2) reaches the FC at most at t^^^^^'' = 

(tVi + 1) i^o + max{7Vi, iVs}] T^; 
the last projection value along 
p2) reaches the FC at most at 



Ni+ N2 + I] uo + Ni 



the paths in 



Overall, the 7r/4 oriented diagonal projections in a quadrant 



are performed after 



iVi + iVs + 1 Uo + max{iVi, 7V2} 



time slots. Again, if the quadrants are processed in a serial 



14 



fashion, the overall number of time slots needed to evaluate 
the diagonal projections sums up to 



TS 



2 [(iVi + N2)iyo + {max{Ni, N2} - 1)] 



2VN [(«! + 02)^0 + maxjai, a2}] — 2 



(11.8) 



Fig. 11 reports, for the sake of clarity, a time diagram summa- 
rizing the timing of projections evaluation when performing 
diagonal projections within a quadrant of the network. 

Then, if we consider a Radon Like scheme comprising P = 
3 projections along the directions '& = 0,1} = Tr/2,'& = 7r/4, 
the overall number of needed transmissions sums up to (cfr. 
(n.3), (II.5) and (II.7)): 



(0,7r/2,7r/4) 



TX 
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2aia2 + 



(ai + a2)max{ai,a2} 



while the number of needed time slots is (cfr. (11.4), (11.6) and 
(II.8)): 



(0,7r/2,7r/4) 



'TS ~ "3^ 

^3 = {oii + 02) (4 + Avo) + max{ai, 02} 

To recap, the Radon Like CS data gathering procedure 
let the fusion center collect all the needed measurement in 
a highly parallelized fashion; far from being optimal, the 
data gathering scheme herein introduced allows a significant 
reduction in both the occupied bandwidth and the consumed 
energy w.r.t. state of the art data gathering scheme as the 
RS introduced in [21]. Further developments of globally 
optimized Radon-like CS data gathering algorithms are still 
under investigation. 
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